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FINITE  ELEMENT  COMPUTER  PROGRAM  TO  ANALYZE 
CRACKED  ORTHOTROPIC  SHEETS 


By 

Chorng-Shin  Chu 
J.  M.  Anderson* 
W.  J.  Batdorf 
J*  A.  Aberson** 


SUMMARY 


The  objective  of  this  study  was  to  develop  a  finite-element  computer  program  that 
performs  a  two-dimensional  elastostatic  analysis  of  plane  anisotropic  homogeneous  sheets 
with  a  through-the-thickness  crack.  The  program  Includes  special  crack-tip  elements  that 
account  for  the  singular  stress  fields  associated  with  crack  opening  (Mode  I)  and  sliding 
(Mode  II)  displacements  at  the  crack  tips.  These  special  crack-tip  elements  provide  a  new 
tool  for  computing  stress-intensity  factors,  and  they  can  be  used  for  predicting  the  crack 
prcpagafion  for  a  damaged  structure. 

Two  types  of  crack  elements  have  been  developed.  One  element  has  8  nodes  and  is 
restricted  to  symmetric  applications;  the  other  element  has  10  nodes  and  is  capable  of  repre¬ 
senting  the  crack-tip  neighborhood  when  both  the  crack  opening  and  sliding  modes  of  defor¬ 
mation  occur.  The  8-node  symmetric  cracked  element  takes  full  advantage  of  the  symmetry, 
so  that  only  one- ha  If  of  a  configuration  which  is  symmetric  about  a  line  containing  the 
crack  need  be  modeled.  However,  the  10-node  unsymmetric  cracked  element  has  much 
wider  usage  in  practical  fracture-mechanics  applications. 

These  cracked  elements  can  exhibit  either  anisotropic  or  isotropic  behavior.  (Ortho- 
tropic  materials  are  considered  a  special  case  of  anisotropic  material.)  A  stress  function  in 
a  half-power  series  form  appropriate  to  the  crack  tip  neighborhood  is  chosen  to  formulate 
the  stiffness  matrix  for  the  anisotropic  cracked  elements,  and  a  Williams*  series  stress  func¬ 
tion  is  used  for  the  isotropic  cracked  elements. 


*Consultant,  Associate  Professor,  School  of  Engineering  Science  and  Mechanics,  Georgia 
Institute  of  Technology. 

**Consuitant,  Assistant  Professor,  School  of  Engineering  Science  and  Mechanics,  Georgia 
Institute  of  Technology. 


The  computer  program  contains  an  axial  element,  a  linear  spring  element,  a  triangular 
element,  and  isotropic  and  anisotropic  cracked  elements.,  Thermahstrain  analysis  capability 
is  included.  Extensive  tests  during  the  development  stage  have  been  made  in  order  to  vali¬ 
date  the  program  against  known  solutions.  To  illustrate  the  capabilities  of  these  two  types 
of  cracked  elements,  a  number  of  selected  sample  problems  are  presented  in  this  report  for 
demonstrations.  The  computer  program  has  proved  to  be  very  efficient. 


INTRODUCTION 


The  finite-element  method  is  one  of  the  most  effective  approaches  available  for  plane 
elasticity  problems  having  irregular  boundaries  or  discontinuous  boundary  conditions.  Early 
efforts  to  bring  this  method  to  bear  on  crack  problems  depended  on  the  use  of  many  conven¬ 
tional  elements  around  the  crack  tip  in  an  attempt  to  represent  the  extreme  stress  gradients 
there.  Stress-intensity  factors  were  estimated  from  results  obtained  with  such  models  either 
by  extrapolating  crack-opening  displacement  (ref.  1)  or  by  numerically  computing  the  vari¬ 
ation  in  strain  energy  with  crack  length  (ref.  2).  Such  conventional  methods  have  been 
found  to  be  limited  and  uneconomical.  For  example,  Oglesby  and  Lomacky  (ref.  3)  have 
indicated  that  the  maximum  permissible  element  size  necessary  to  ensure  acceptable  accu¬ 
racy  (5%  error  or  less)  is  on  the  order  of  1/500  of  the  crack  half  length. 

To  circumvent  this  economic  problem,  development  and  research  efforts  have  turned 
toward  formulating  elements  which  contain  the  crack-tip  stress  singularity.  These  special 
singularity  elements,  usually  referred  to  as  cracked  finite  elements,  represent  a  significant 
improvement  in  both  accuracy  and  economy  In  comparison  with  conventional  methods.  Many 
cracked  finite  elements  developed  to  date  (ref.  4-8)  Incorporate  only  the  singular  term  in 
the  series  expansion  for  the  crack-tip  stress  field.  Admittedly,  this  term  dominates  all 
others  near  the  crack  tip,  but  to  guarantee  that  the  nonsingular  contributions  are  comparably 
negligible,  the  neighborhood  represented  by  the  cracked  element  must  be  quite  small,  and 
the  problem  of  economy  arises  again.  Moreover,  the  neighboring  conventional  elements  in 
such  instances  are  drawn  very  near  the  crack  tip  again,  and  concern  over  their  capability 
to  represent  the  stress  field  adequately  has  led  some  investigators  to  introduce  special 
"border  elements"  having  a  higher  degree  of  sophistication  than  that  routinely  required  for 
plane  elasticity  problems, 

Wilson  (ref.  9)  has  developed  a  cracked  finite  element  that  makes  use  of  the  first  four 
terms  in  the  expansion  for  the  crack-tip  stress  field.  He  reports  accurate  results  when  this 
element  is  used  in  conjunction  with  a  fairly  modest  number  of  conventional  elements. 
Wilson*s  element,  however,  has  the  disadvantage  of  being  semi-circular  and  hence  is  some¬ 
what  awkward  to  use  with  conventional  elements,  which  almost  always  have  straight  bound¬ 
aries.  Moreover,  Wilson*s  element  (as  well  as  some  others  previously  referenced)  has  fewer 
degrees  of  freedom  than  are  needed  for  independence  of  the  nodal  displacements.  This 

requires  that  the  stiffness  matrix  of  the  cracked  element  receive  special  attention  in  form¬ 
ing  the  stiffness  matrix  of  the  assembly. 
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Early  In  1973,  the  Lockheed-Georgla  Company  completed  the  development  of  isotropic 
versions  of  two  high-order  cracked  finite  elements,  i.e.,  elements  that  incorporate  many 
terms  in  the  expansion  for  the  crack-tip  stress  field.  This  feature  permits  very  accurate  esti¬ 
mates  of  stress-intensity  factors  with  relatively  coarse  finite-element  grids.  Both  high-order 
cracked  elements,  one  for  symmetric  (Mode  I)  applications  and  one  for  unsymmetric  (Mode  I 
and  Mode  II)  applications,  have  a  perfect  balance  between  actual  degrees  of  freedom  and 
the  number  of  nodal  displacement  components.  Thus,  the  numerical  analyst  adds  the  cracked 
element  to  an  assembly  in  exactly  the  same  way  that  he  adds  a  conventional  element.  The 
shape  of  each  element  was  chosen  to  make  it  fit  conveniently  in  models  making  use  of  the 
most  widely  used  membrane  element:  the  constant-strain  triangle.  The  remarkable  accuracy 
obtained  with  both  elements  in  isotropic  applications  using  very  coarse  finite-element  repre¬ 
sentations  (ref.  10)  would  seem  to  indicate  that  the  periodic  concern  of  some  investigators 
over  the  displacement  incompatibility  that  exists  at  the  interface  of  the  high-order  cracked 
element  and  a  conventional  element  is  largely  academic. 

The  continually  increasing  use  of  fiber-reinforced  composites  in  aerospace  applications, 
coupled  with  a  growing  confidence  in  the  ability  of  linear  elastic  fracture  mechcriics  (LEFM) 
to  predict  the  growth  rate  and  stability  of  cracks  in  isotropic  materials,  has  lately  resulted  in 
considerable  interest  in  the  prospects  of  successfully  applying  LEFM  to  anisotropic  materials 
(ref.  11).  This  report  contains  the  analytical  background  and  an  account  of  several  numeri¬ 
cal  modifications  and  additions  required  to  effect  the  following  principal  extensions  of 
Lockheed-Georgia's  existing  analysis  capability  for  cracked  isotropic  structures: 

o  Axial  element 

o  Anisotropic  constant-strain  triangle 
o  Anisotropic  cracked  elements 

o  Automatic  node  sequencing  to  minimize  bandwidth 
o  Input-data  generator 
o  Thermal  stress  analysis  capability 


SYMBOLS 


ipf 

{p} 

[k] 

hi 


Vector  of  forces  for  uncoupled  structure 

Vector  of  displacements  for  uncoupled  structure 

Vector  of  thermal  forces  for  complete  restraint  of  uncoupled  structure 

Uncoupled  stiffness  matrix 

Vector  of  displacements  for  coupled  structure 
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(Q  1 

[*] 

[K] 

IM 

[U] 

)“( 

X/  y,  z 
r,  0 


Vector  for  forces  for  coupled  structure 
Compatibility  matrix  for  displacements 
Coupled  stiffness  matrix 

Vector  of  thermal  forces  for  complete  restraint  of  coupled  structure 
Stress  vector 
Strain  vector 

Elastic  matrix  for  Hooke's  law  in  local  axis 
Elastic  matrix  for  Hookes  law  in  material  axis 
Transformation  matrix  for  elastic  coefficients 
Vector  representing  thermal  expansion  coefficients 
Rectangular  coordinates 
Polar  coordinates 
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T 


Strain  energy 
Temperatu  re 


Stress  function 
Half  crack  length 

Total  height  of  cracked  tension  plate 
Total  width  of  cracked  tension  plate 

=  W/2 

Crack  length  aspect  ratio 


BASIC  EQUATIONS  FOR  FINITE  ELEMENT  PROGRAM 


Matrix  Equations  for  Stiffness  Method 

The  finite  element  analysis  program  is  based  on  the  direct  stiffness  method.  With  this 
approach,  a  complex  structure  is  idealized  as  an  assembly  of  simple  elements;  for  example: 
axial  elements,  triangular  membrane  elements  (isotropic  and  anisotropic),  and  8-  and  10- 
node  cracked  elements.  Each  of  these  elements,  separately,  can  be  analyzed  without  diffi¬ 
culty.  By  expressing  the  various  mathematical  relationships  in  matrix  form,  it  is  possible  to 
assemble  automatically  a  large  number  of  discrete  elements  to  simulate  almost  any  complex 
structure.  The  basic  matrix  operations  that  are  executed  by  the  program  are  illustrated  here. 
These  equations  include  thermal  effects. 

In  matrix  form,  the  uncoupled  force-displacement  equation  for  a  single  element  or  the 
entire  structure  can  be  written  as 

{p}=[k]{p}+{p^j  (,) 

where  P  and  p  are  vectors  of  forces  and  displacements,  respectively,  and  k  is  the  uncoupled 
stiffness  matrix.  The  term  is  a  vector  of  thermal  forces  equivalent  to  the  forces  produced 
when  each  element  Is  completely  restrained. 

From  compatibility  considerations  we  can  define  il)  such  that 


where  q  is  a  vector  of  node  point  displacements  for  the  coupled  structure.  In  other  words 
the  p  vector  is  referred  to  the  local  uncoupled  system  and  q  to  the  coupled  structure  de"" 
fined  in  some  global  reference  frame.  The  matrix  0  consists  of  direction  cosines  and  is 
determined  by  the  topology  of  the  structural  model.  The  vector  q  can  then  be  expressed  in 
terms  of  unknown  displacements  q^  and  known  displacements  From  the  principal  of 

virtual  work  it  can  be  shown  that 


(3) 

{q}=  [K]{q}+|Q„} 

(4) 

Equation  (4)  is  then  the  force^displacement  equation  for  the  coupled  or  global  system, 
where  Q  and  q  are  forces  and  displacements,  respectively.  The  vector  Q^|.  represents 
forces  for  complete  restraint  for  the  global  system.  Equation  (4)  can  be  expanded  as 


where  Qa  and  qg  are  applied  forces  and  unknown  displacements;  and  Qj^  and  q^^  are 
reaction  forces  and  known  displacements.  To  determine  the  elastic  displacements  for  a  hot 
structure  it  is  then  necessary  to  solve  the  following  set  of  linear  equations: 


Once  the  displacements  are  known,  the  loads  in  each  discrete  element  can  then  be  deter- 
mined  from  the  following  equation 


The  reaction  forces,  along  with  the  equilibrium  check,  can  be  determined  from  Equation  (8) 


References  (12)  and  (13)  discuss  the  various  concepts  of  matrix  analysis  of  structures. 
The  derivations  of  the  element  stiffness  matrices,  except  for  the  cracked  elements,  are  out¬ 
lined  in  the  Appendix. 


Transformation  of  Elastic  and  Thermal  Coefficients 

For  a  two-dimensional  elastic  plate  the  relationship  between  stress  and  strain  can  be 
written  as 


{''}=  k]{'} 


For  an  isotropic  material  and  plane-stress  conditions, 

E  vE  ^ 
r — JT  - s-  0  1 


E 

vE 

l-i/^ 

vE 

E 

1-y^ 

]-u^ 

In  the  general  case  of  anisotropic  materials  the  only  requirement  is  that  the  matrix  of  co- 
efficients  be  symmetric,  that  is 
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AAA 
11  12  16 


K‘12  ^22  ^26 


AAA 
16  ^26  ^66 


The  six  coefficients  are  defined  with  respect  to  a  material  axis  whose  orientation 
usually  does  not  correspond  to  the  principal  axis  of  the  structural  element.  The  transforma¬ 
tion  matrix  used  within  the  program  to  rotate  the  elastic  coefficients  was  expressed  in  terms 
of  direction  cosines.  It  is  very  important  to  note  that  all  the  angles  are  measured  from  the 
local  axis  to  the  material  axis  (see  figure  1). 

The  transformation  for  the  elastic  coefficients  is  then 


=  U  ■  A 


where 


2 

cos  a 


U  =  cos  /3 


2/? 
cos  p 


2 

cos  0! 


-2cos0i  cos^  2cosQ:  cos 


cos  a  cos/S 


~CO$Cl  cosp 

2  2 
cos  0!  “  cos  /S 


and  for  the  thermal  expansion  coefficient  vector 


In  the  computer  program  all  these  transformations  are  performed  automatically  for  each 
element.  It  is  only  necessary  to  indicate  the  direction  of  the  material  axis  by  specifying 
two  points  on  the  axis.  Several  material  axes  can  be  specified  for  a  single  structural  model 
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TWO-DIMENSIONAL  ELASTIC  SOLUTIONS  FOR 
CRACK-TIP  STRESSES  AND  DISPLACEMENTS 


The  field  equations  of  elasto  statics  appropriate  to  isothermal  deformation  under  plane- 

stress  conditions  (a  =  T  =  T  =0)  with  zero  body  forces  are  given  below  and  represent  a 
z  yz  zx 

basis  for  all  the  equations  developed  in  this  section. 


Bu  9u  bu  bu 

e  =  — ^  ,  e  =  and  Y  =  -r-”  + 

X  3x  y  by  xy  by  bx 

2  2  2 
be  b  e  by 

_ 1  = 

^.2  ^2  bxby 


°n  °]2  °]6 


°]2  °22  °26 


a, ,  a 


16  “26  “66 J  \  xy 


ba  bT 

bx  by 


ba  bT 

_X+_i^  =  0 

by  bx 

For  plane-strain  applications  (e  =Y  =  Y  =0),  the  constants  a.,  in  (16)  must  be 

'  z  yz  zx  1 1 

replaced  by  their  plane-strain  counterparts  bj.  given  by 

(i,i  =  1,2,6) 

1]  i|  0^3 


The  form  of  the  equilibrium  equations  (17)  implies  the  existence  of  a  stress  function 
U(x,y)  such  that 

2  2  2 
„  U  ^  _b^U  .  ,  _  b^u 

3/2  y  3x2  xy  3x3y 
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By  using  (19)  to  eliminate  the  stress  components  in  (16)  and  substituting  the  resulting  strain 
components  into  the  compatibility  condition  (15),  we  find  that  the  stress  function  must  satisfy 

4  4  4  4  4 

9  U  ^  5  U  ^  B  U  ^  B  U  ,  b  U  _ 

°22  .  4  ^°26  3  ^^°12  ^66^  2^  2  "  ^°16 ,  ^  3  °1K  4  ° 

OX  ox  oy  3x  oy  9xdy  9y 

Equation  (20)  may  be  written  as  the  product  of  four  first-order  operators 

'^i""^"^iSr  (21) 

If  \Jl^,  (ig  and  are  roots  of  the  characteristic  polynomial, 

°1 1^"^  "  ^'=’16^^  °66^^  ■  2°26^  °22  "  ^ 
then  (20)  may  be  written  as 

D^D2D3D^U(x,y)  =0  (23) 

The  roots  of  (22)  occur  in  complex-conjugate  pairs  and  can  be  shown  by  energy  considerations 
always  to  have  non-zero  imaginary  parts.  The  designation  u^,  ,l^,  and  shall  henceforth 

be  used  for  the  roots  of  (22),  and  it  will  be  taken  for  granted  that  IJ.  -  |J..  0.  The  general 

solution  of  (23)  for  U(x,y)  can  now  be  written  by  repeated  integration,  but  the  form  it  takes 
depends  on  whether  or  not  the  roots  of  (22)  are  distinct.  A  discussion  of  cases  follows. 


Distinct  Roots  (n  ^  ^  IJ.2) 

For  distinct  roots,  the  general  solution  of  (23)  for  real  U(x,y)  is 

U(x,y)=2Re[U^(zp  +  U2(z2)] 

where  and  U2  are  arbitrary  functions  of  the  complex  variables 
=  X  +  pt^y  and  Z2  =  x  +  i-i2y 

respectively.  Upon  substituting  (24)  into  (19)  the  following  expressions  are  found  for  the 
stress  components. 

=2Re[^i^U^'(zp+^i2^J..^^^)J 


a  =2Re[Uf(zp+U''(z2)] 
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and  =  -2Re[ti^  U J(zp  +  U2U2(z2)l 

in  which  prime  (')  denotes  differentiation  with  respect  to  the  parenthetical  argument.  It  may 
be  verified  by  differentiation  that  the  displacement  components  corresponding  to  (26)  are 
g  iven  by 


u^=2Re[p,U;(zp+P2U-(z2)] 

and  u  =2Re[q^U'^(zp +q2U2(z2)] 


provided 


’’l  '^°I2  '  '’2  '^°12  ■  “id't 


and  qj=a,2M,  +  — -=2^ 


'  ^2  12  2  ^2  26 


and  U2  components  of  the  stress  function  are  taken  to  be  half-power  series,  i.e,, 

00  n-2  “  n-2 

Ui'f-,)  =  2  z,2  and  U''(z  )  =  ^  {29] 


These  will  now  be  shown  to  satisfy  boundary  conditions  appropriate  to  free  crack  faces  by 
properly  relating  the  complex  coefficients  C  and  D  .  From  (26),  the  stresses  corresponding 
to  a  typical  term  in  (29)  are 

n-2  n-2 ' 

2  T~  2  ~T' 

G  =2Re  ii,C  z,  +H^  D  z,^ 

X  .  I  n  1  2  n  z  J 


G  =  2  Re  C  z ,  +  D  zJ 

y  L  n  1  n  2 


and  T  = -2Re  t-i,  C  z,  D  z. 

xy  L  1  n  I  2  n  2 


Upon  considering  the  crack-tip  neighborhood  and  coordinate  system  shown  In  figure  2, 
we  can  write 


G  (x,0)  =  T  (x,0)  =  0  for  X  <0 
y  xy 
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as  the  crack-face  boundary  conditions.  From  (30),  these  require 
Re 


and 


Im  ^^n 


Re 


(C..D„)=0  ,whennh-='’ 


(32) 


Equations  (32)  are  satisfied  if  we  set 


and 


(33) 


in  which  i  is  the  imaginary  unit  and  and  are  arbitrary  real  constants.  Using  (33)  to 
write  C  and  D  in  terms  of  A„  and  B„,  we  obtain 


and 


^  ,.,n+1 

=n  <'>  irr^n- 


, ,  B  -  fj.  A 

D  _ Ln 

n  ^2  ” 


(34) 


From  (26),  (27),  (29),  and  (34),  we  are  now  able  to  write  the  crack-tip  stresses  and  dis¬ 
placements  corresponding  to  a  typical  term  in  the  and  U2  series. 


a  =  2Re 

X 


(35) 


and  T  =  -2Re 


(36) 
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Repeated  Roots  (M-^  -\^2  “ 

In  this  cose  the  operators  in  (23)  are  not  distinct;  each  is  once  repeated;  i,e.. 


where 


Dj[£u(x,y)  =  0 


D  and 

H  Sy  3x  U  3y  ox 


Repeated  integration  now  gives 

U(x,y)=2Re[U^(z^)  +I^U2(zi)] 
as  the  general  real  solution  of  (37).  In  (39), 


z  =x+Hy  and  z^=x+M-y 


Expressions  for  the  stress  components  corresponding  to  (39)  are  found  from  (19)  to  be: 

!7^  =2Re  [iJ  U;'(z^)  +2^fIU^(zp  z^U^'(zj)] 

=2Re[UJ(zp +2U^(z^) +I,U^'(z^)]  (41) 

and  =-2Re[|aUJ(zp  +  (^  +il)  U2(Z|)  +Iiz^U2(z^)] 

Moreover,  the  displacements, 

u^=2Re[p^U'(z^)  +P2U2(z,)  +P3li4(zp] 

and  u^  =2Re[q^U’ (zp  +q2U2(z,)  +q3l,U^(zp]  (42) 

can  be  shown*  to  be  consistent  with  the  stresses  (41),  the  stress-strain  equations  (16)  and  the 
strain-displacement  equations  (14)  provided 

*ln  verifying  (^)  subject  to  (14),  (16),  (41),  and  (43)  certain  identities  must  be  used  appro¬ 
priate  to  Hand  \i  as  repeated  roots  of  (22).  These  are: 

^°12  ^°66  “  °11^^  +4|JH+ii) 

and  a22  “ 
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Pi  P3  ^  °11  ■^"’12  ~^°]6 


P2  -  M-(2m.  -  m)  ^  +  a^2  “  1-^ 


Pi  P3  °12^  ^  H  “  °26 


and  q  =a  ^ +  — (2  --)- a-, 

2  12  (i  \  fi/  26 


Again,  to  develop  a  half-power  series  analogous  to  (29),  we  take 


=  I  A 


in  which  Cp  and  are  complex  constants  and  n  is  anticipated  to  be  a  real  integer. 
Following  the  same  procedure  used  for  the  distinct-root  case,  we  substitute  (44)  Into  (41)  and 
consider  crack-face  boundary  conditions  (31).  This  leads  to  the  requirement  that  C  and  D 
ivplcal  !y  SQi  isfy  n  n 


/-  _.n+l^ 

C  +— —  D  =  I  A 
n  2  n  r 


and  fiC  +(li  D  =  B 

n  \  2  /  n  f 


In  which  A  and  B  are  real  constants.  The  simul; 

.  I  I  n  n 
yields 


aneous  solution  of  (45)  for  C  and  D 

n  n 


^n  — 


and  D  = - (B  -  laA  ) 

n  -  n  n' 

pi  -  la 


When  these  expressions  are  incorporated  in  (41)  and  (42),  the  following  formulas  for 
stresses  and  displacements  evolve  for  a  typical  term  of  the  series: 
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f  n-2  r  /  - 

.n+1  -sr  o  I  _ 

o  =  Re  - z,  A  M.  (n|a- 2ia  -  |i(n  -  2)  — 

X  —  I  n  \  z, 

u  -n  \  1 


+  B^li  I4|i  -  (i(n  +2)  +|a(n  “  2)  - 


•n+l 

Q  =Re  Z - {ia(n  -  4) +2m.  -  (i(n  -  2)  —  j. 

^  \  1/ 


I  i"*'  2  /  -  '‘l 

T^y=Re  —  z,  u  A_,(n  -  2>(  -  -  1 1  +  U-2^-^(n-2)- 


22  r  /  - 

\  =  Re  - - z^  -  lp^(2u  +n^)  -  2p24  - 

I  \  1 


+  B  (2p2-p,(n+2)+np,- 


I2  r  /  - 

^  H-|i  L  \  1 


+  B^(  2q2  -  q^(n  +2)  +nq^ 


The  repeated-root  cose  includes  the  stress  and  displacement  functions  appropriate  to  the 
crack-tip  neighborhood  in  an  isotropic  material.  For  an  isotropic  material  and  plane  stress 
conditions. 
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°11  °22  E 


a 


12 


V 

E 


"16  °  "26  '  “ 


and 


_2(1  +v) 


66 


With  these  simplifications,  the  characteristic  polynomial  (22)  reduces  to 

4  2 

li  +2n  +1=0 

which  has  as  a  solution  the  repeated  roots 
U  =  ±  i 
Consequently, 

U(x,y)  =2Re  [U^(z)  +  zU^{z)] 

in  which 

z  =  X  +  iy  and  z  =  x  -  iy 
Upon  integrating  (44)^,  we  find 

n+2  ^ 

U,(z)  =  C*  z  ^  and  U_(z)  =  D*  z^ 

1  n  2  n 

as  components  of  the  isotropic  form  of  U(x,y).  Combining  (52)  and  (54),  we  find 

n+2  n+2  n  n 

U(x,y)=C*z^  +^z  ^  +z  D*  z^  +z  D*  z^ 
n  n  n  n 

In  polar  form  (see  figure  2), 


(49) 


(50) 


(51) 


(52) 


(53) 


(54) 


(55) 


In  (54)  it  Is  convenient  to  use  C*  and  D*  to  represent  4C  /n(n+2)  and  2D  /n,  respectively. 
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n  n  .n£  n 

2  2*2  2  /  ne  n9\ 

z  =r  e  =r  (  cos  ^  +  i  sm-^j 

n  n  .n9  n 

-2  2  “'2  2  /  n9  .  .  n9 

z  =r  e  =r  \^cos — 1  sin 


so  that  (55)  becomes 


U(r, 6)=C*r^  ^cos^j  +  1^  9+isin^j  +  1^0 


+  C*  r  {  cos 
n  \ 


j  +  1  j  9  -  i  sin(|  +  9 


+  D*  r  ^  f  cos  ( J  "  ^)  9  +  i  sin  ^  j  9 


+  D*  r  ^  ^  1)  9 


=  2r  ^  [Re(C*)cos(j  +  l)  9  -  lm(C*)  sin  (  j  +  l)  9 
+  Re(D*)  cos^j-  9  -  lm(D*)  sin^j  -  1^9^ 


which  is  a  typical  term  in  the  familiar  Williams  series  of  stress  functions  appropriate  to  the 
crack-tip  neighborhood  in  an  isotropic  material  (ref.  14,  15).  Imposition  of  the  crack-face 
boundary  conditions  (31)  leads  to 


Re(D*)  = 
n 


±1 

j + (-1)' 


Re(C*) 


lm(D*)  = 
n 


— = - lm(C*) 
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Upon  substituting  (58)  into  (57),  we  find  the  following  expressions  for  stresses  and  displace¬ 
ments  in  polar  form: 


and 


,  1  8U  1  a  U 

^r  r  3r  2  ^2 

r  30 


n-2 

=  lRe(C*)r2  n 


-  (n  +  2)  cos  ( — :r— — — — — —  (n  -  6)  cos  i )  9 


n-2 

+  l|m(C*)r  ^  n 
2  n 


n  +2(-1) 


-  (n  +  2)  sin  f— )  9  +  — - (n  -  6)  sin  -  1 0 


n-2(-iy 


a^u 


ar 


n-2 

=  iRe(C*)  r  2  n 
2  n 


!  J.n\  n+2  /n-'> 

(n  +2)  cost  — — ) - cos 


n-2 

+  ;7  lm(C*)  r  ^  n(n  +2) 
Z  n 


a  /I  au 


^  ^  n  +2(-l)" 


./'n+2\Q  n+2  ,/n-2,Q 

sin  ^ — 1  9 - sin  I  — r: —  1  9 


n  -  2(-l) 


r9  Srlr  S9 


n-2 

=  lRe(C*)r  ^  n 
2  n 


(n  +2)  sin  f— )  9 - HJl? —  (n  -  2)  sin  )  9 


n  +2(-l)' 


n-2  ^ 


+  ^  lm(C*)  r  ^  n 
2  n 


-  (n  +  2)  cos  ( —5—)  9  +  — —  (n  -  2)  cos  i  -  ^  1  9 

^  n-2(-1)" 


18 


u  =i^^Re(C*)r^ 


(n  +  2)  cos  (4^)  9  -  (6  --^  .  „)  cos  (l^i)  9 

n+2(-l) 


1  +v 
E 


|m(C*)r^  (n +2)sin(^^^e 


.  n+2  (,  8v  \  .  /n-2\a 


1  +v 


E 


R.  (C*)  [(n  «)  s!n  (ll2)  9-  +  ")  (=^)  « 

n+2(-l)  J 


+  i^  |m(C*)r^  .(„+2)cos(!l|2)e+— 2-^(6-^+n 
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FORMULATION  OF  CRACKED  ELEMENT  STIFFNESS  MATRICES 


Two  types  of  cracked  elements,  as  shown  in  figure  3,  have  been  developed  and 
implemented,  because  many  fracture  mechanics  problems  are  symmetric  about  the  plane  of 
the  crack.  One  formulation  takes  only  the  symmetric  terms  in  the  series  solution  and, 
hence,  is  applicable  only  to  symmetric  problems  (K||  =  0);  the  other  formulation  makes  use 
of  both  symmetric  and  anti-symmetric  terms  and  is  applicable  to  unsymmetric  or  mixed-mode 
problems  (K|  and  K||). 

The  coordinate  system  of  an  8-node  symmetric  element  has  its  origin  at  the  crack  tipo 
It  Is  rectangular  in  shape  with  a  three-to-one  aspect  ratio.  Placement  of  the  nodes  relative 
to  the  rectangle  is  pre-determined  with  a  node  at  each  comer  plus  nodes  at  the  one-third 
points  of  each  of  the  long  sides.  The  lower  side  (nodes  6,  7,  8,  and  1)  is  coincident  with 
the  crack  direction  and  presumed  axis  of  symmetry.  Nodes  6  and  7  are  on  the  free  crack 
face.  Nodes  8  and  1  are  on  the  prolongation  of  the  crack.  They  are  constrained  rigidly  as 
to  vertical  displacement  and  are  free  of  shear  forces  -  conditions  that  are  consistent  with 
symmetry. 

The  8-node  symmetric  element  has  16  displacement  degrees  of  freedom,  two  per  node 
corresponding  to  the  in-plane  displacement  components.  Thus,  it  incorporates  the  first  13 
symmetric  terms  of  the  series  associated  with  rigid-body  morion  in  the  plane.  These  13 
coefficients  and  the  three  rigid-body  parameters  are  referred  to  as  the  16  generalized  co¬ 
ordinates  of  the  cracked  element.  The  stresses  and  displacements  corresponding  to  these  16 
generalized  coordinates  are  evaluated  on  the  boundary  of  the  element.  Products  of  stress 
and  displacement  contributing  to  boundary  work  are  performed  and  integrated.  The  result  is 
a  homogenous  quadratic  form  in  the  generalized  coordinates,  and  the  coefficient  of  each 
term  is  an  element  of  the  cracked  element  stiffness  matrix  with  respect  to  the  generalized 
coordinates.  Once  the  stiffness  matrix  with  respect  to  generalized  coordinates  is  deter¬ 
mined,  the  stiffness  matrix  with  respect  to  nodal  displacements  is  formed  using  the  series 
(with  rigid-body  terms)  to  write  nodal  displacements  in  terms  of  the  generalized  coordinates. 

The  lO-node  unsymmetric  cracked  element  as  shown  in  figure  3  is  square  with  equally 
spaced  nodes  around  its  boundary.  Like  the  symmetric  cracked  element,  the  shape  of  the 
element  and  relative  location  of  the  nodes  were  chosen  to  provide  modeling  convenience. 
The  generalized  coordinates  correspond  to  the  first  9  symmetric  terms  and  first  8  anti-sym- 
metric  terms  of  the  series,  plus  the  3  rigid-body  displacement  parameters.  The  stiffness 
matrix  was  again  generated  by  integration  around  the  boundary. 

To  formulate  the  cracked-element  stiffness  matrix,  it  Is  convenient  to  introduce 
dimensionless  variables  for  both  anisotropic  and  isotropic  cases. 
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Anisotropic  Case 


The  dimensionless  variables  for  the  anisotropic  case  are  given  as 


z  .j  /  A  / 

=  z^/ A 

P,  =P,/a,|  , 

Ps'PsAll 

Q,  , 

Q2  =  q2/'>,| 

C  =  a,  ^  A  A  , 

nil  n 

5-1 

c 

CO 

<3 

o 

1! 

c 

Q 

S  Q-i  ,  CT  , 

X  1  1  X 

S  a . .  ^  , 

y  ' '  y 

U  =u  /a 

X  x" 

U  =  u  /A 

y  y 

(61) 


where  A  is  a  characteristic  length  and  is  taken  to  be  the  distance  between  nodes  on  the 
cracked  element.  In  terms  of  the  dimensionless  variables  in  equation  (61),  equations  (35) 
and  (36)  take  the  form  below 
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and 


A  (  /-xri  +  1 

i  Re  U;^ - 

n 


Q 


Equations  (62)  and  (63)  are  the  basic  equations  needed  to  form  the  stiffness  matrix  of 
the  anisotropic  cracked  element.  The  strain  energy  V  stored  In  the  cracked  element  can 
be  computed  by  numerically  integrating  the  work  of  the  surface  tractions  around  the  bound' 
ary.  This  yields 


_  h  A  T , 

2V  -  -  q  kq 

a, ,  ^  ^ 


in  which  h  =  the  uniform  thickness  of  the  cracked  element 
q  -  the  column  matrix  of  generalized  coordinates 
q^  =  the  transpose  of  q 

k  —  the  stiffness  matrix  with  respect  to  generalized  coordinates 

Let  D  denote  the  column  of  dimensionless  nodal  displacement  components  in  the 
elemental  cartesian  coordinate  system.  The  matrix  C  in 

D  =  Cq  ( 

is  obtained  by  evaluating  equations  (63)  at  the  nodes,,  The  inverse  of  C  may  be  found 
numerically  with  the  result 

q  =  C  ( 

Thus,  from  (64), 

2V  =  D^(C“^  —  k  C“^D  ( 

°n  ^ 


or  in  terms  of  D*,  the  column  of  dimensional  nodal  displacements 
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2V  =  D*^(C'^  —  kC’^  D* 

Q  <«  T 


K  =  —  (C  )  k  C 
°11 


Once  the  stiffness  matrix  of  the  cracked  element  is  formed,  its  incorporation  into  the 
stiffness  matrix  of  an  assembly  follows  exactly  the  procedure  used  for  conventional  elements. 
The  same  approach  was  followed  for  the  case  of  repeated  roots. 


Isotropic  Case 

The  dimensionless  variables  introduced  for  an  isotropic  case  are  given  below: 


A  ' 


G  ' 


n  2G 


R  A 

S  =  ^ 
9  G 


Re  (C*) 
n 


s 

Re  G 


5“’ 

where  A  is  also  a  characteristic  length  as  defined  before. 

In  terms  of  the  dimensionless  variables  in  equation  (70),  equations  (59)  and  (60)  take 
the  following  form: 

I  "  ’ 

S|^(R,  9)=  nR  |s*  j^- (n +2)  cos  (^+ 1)  9  +  f  (n)  (n  -  6)  cos  (^-  l)ej 

+  a*  j^g(n)(n  +  2)  sin  (^  +  i)  9  -  (n  -  6)  sin  (^  -  1)  9  | 
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Sg  (R,  9)  = 


n  =  1 


2  ^  r  1 

n  (n  +  2)  R  Is*  cos  (^  +  1)0-  f(n)  cos  (^  -  1)  0 


+  a  *  ”  g{n)  sin  (:^  +  ])  0  +  sin  (tt  -  1)  0 
n  ^  z 


SRg(R,e)  = 


n  =  1 


1  v>} 

n  R  l^n  (^+1)9"  f(n)(n  -  2)  sin  -  1) 


+  a*  g(n)(n  +  2)  cos  (^  +  1)  0  -  (n  -  2)  cos  (^  -  1)  0 


Up(R,  0)  -  R  ) 

^  n=  ]  ^ 


®  £  (  r 

^  j^2  I  s*  -  (n  +  2)  cos  g  +  1)  0  -  f(n)(6-8  ?-n)  cos  (J  -  1)  0 


a*  (n +2)  g(n)  sin  (J+ 1)  0 +(6-8  4-n)  sin  (J-  1)  0  | 


00  j 

,(R/  6)  =  ^  ^  R  I  s*  (n  +  2)  sin 

n  =  1  (  "  L  ^ 


^  +  1)0“  f(n)  (6=84  +  n)  sin  (^  ~  1)  0 


a*|^(n  +  2)  g(n)  cos  (^  +  1)0-  (6-8^  +  n)  cos  (^  -  1)0  | 


in  which 


f(n)  = 


5  +  (-1)" 


g(n)  = 


i>  for  plane  strain 

V  for  plane  stress 

l+l^ 
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Equations  (71)  and  (72)  are  the  basic  equations  needed  to  form  the  stiffness  matrix  of  the 
isotropic  cracked  element, 

it  is  noticed  that  equations  (72)  for  the  displacement  components  have  been  v/ritten  so 
as  to  keep  the  parts  independent  of  ^  and  distinct  from  the  parts  dependent  on  4*  This 
distinction  was  made  to  permit  the  storage  of  the  stiffness  matrix  as  the  sum  of  two  matrices 
that  are  each  internally  independent  of  material  parameters. 

By  numerical  integration,  the  strain  energy  V  stored  in  the  cracked  element  can  be 
written  in  terms  of  the  stiffness  matrices  k-j  and  k2,  i.e., 

2V  =  GhA^q^  (k^  +Sk2)  q  ^^4) 

in  which  h  is  the  uniform  thickness  of  the  element.  Following  the  same  approach  as  for 
the  anisotropic  case,  equation  (74)  can  also  be  written  in  terms  of  D*,  the  column  matrix 
of  dimensional  nodal  displacements, 

T  1  T  , 

2V  =  D* '  (C" ')  Gh  (k^  +  C  D*  (75) 


So  that 


K  =  Gh  (C 


T 


-1 


(76) 


Each  of  the  stiffness  matrices  of  these  isotropic  cracked  elements  is  stored  as  the  sum 
of  two  matrices  that  are  independent  of  the  size  and  material  properties  of  the  cracked 
element.  This  means  that  the  boundary  integration  mentioned  previously  does  not  have  to 
be  repeated  for  each  application.  Consequently,  the  efficiency  of  the  analysis  program  in 
conventional  finite“element  applications  is  not  diminished  at  all  by  the  addition  of  the 
cracked  element.  The  displacement  incompatibility  that  exists  between  nodes  where  the 
cracked  element  interfaces  with  a  conventional  element  seems  inconsequential  in  light  of 
the  exceptional  accuracy  obtained  in  many  and  varied  applications,  including  the  sample 
prob  lems . 
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STRESS-INTENSITY  FACTORS  AND  STRAIN-ENERGY  RELEASE  RATES 


Stress- Intensity  Factors 

- 1  /2 

The  leading  terms  in  equations  (35)  and  (59)  contain  the  singularity  r  '  ;  all  subse¬ 
quent  terms  are  non-singular.  The  coefficients  A-j  and  B]  (C|*  for  isotropic  case)  are 
related  to  the  opening  and  sliding  mode  stress- intensity  factors  Ki  and  Kii  by  the  following 
formulas:  " 


Anisotropic  Case: 

K,  =  a^(x,  o)  =  -2/^ 

(77) 

K|i  '^xy(xr  o)  =  2/^ 

Isotropic  Case: 


K,  =  C7g  (r,  o)  =  6  Re  (C^*) 

K||  =  y  2ttr  “^rOir,  o)  =  -2/  27r  Im  (C^*) 


(78) 


Strain-Energy  Release  Rates 

The  strain-energy  release  rates  are  related  to  the  stress- intensity  factors  by  the  elastic 
coefficients  in  the  following  fashion: 

Gj=  cK.^  i  =  |,  II  (79) 

where  "i"  indicates  the  mode  number.  The  total  strain-energy  release  rate  is  the  summation 
of  energy  rate  at  each  mode  based  upon  linear  superposition.  The  elastic  coefficients  c 
that  relate  energy  rates  to  stress-intensity  factors  are  listed  in  Table  I. 
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ELASTIC  COEFFICIENTS  RELATING  ENERGY  RATES  TO  STRESS- INTENSITY  FACTORS 


BRIEF  DESCRIPTION  OF  CRACKED  ELEMENT  COMPUTER  PROGRAM 


The  computer  program  is  an  a  I  h FORTRAN  program  which  uses  only  conventional  l/O 
(FORTRAN  READ  and  WRITE  statements)  and  the  standard  library  routines,  it  contains  the 
following  finite  elements: 

•  Axial  element 

•  Linear  spring  element 

•  Plane-stress/strain,  isotropic  and  anisotropic,  triangular  element 

•  Eight-node  symmetric  isotropic  cracked  element 

•  Eight-node  symmetric  anisotropic  cracked  element 

•  Ten-node  unsymmetric  Isotropic  cracked  element 

»  Ten-node  unsymmetric  anisotropic  cracked  element 

The  computer  program  consists  of  a  main  program  called  CRAK,  which  allocates  the 
core  and  calls  sequentially  four  main  driver  subroutines  (INPUT,  KFORM,  CHOL, 
OUTPUT),  and  32  subroutines.  These  subroutines  have  been  overlayed  to  reduce  the  core 
required  for  program  code,  permanent  data,  and  program  table.  The  segment  structure  of 
the  program  Is  shown  in  figure  4.  The  name  and  function  of  each  subroutine  are  given  in 
Table  II. 

The  input  segment  of  the  program  contains  data  generator  features  to  minimize  the 
labor  required  to  Input  a  model.  During  the  input  sequence,  a  banding  algorithm  Is  used 
to  resequence  the  node  numbers  to  minimize  the  band-width  of  the  structural  stiffness 
matrix  which  is  assembled  in  the  KFORM  segment.  This  banding  algorithm  is  quite  effec¬ 
tive  and  helps  to  minimize  execution  times.  A  banded  Cholesky  decomposition  procedure 
(segment  CHOL)  is  used  to  solve  the  system  of  equations.  The  computer  program  organiza¬ 
tion  is  especially  characterized  by  its  flexibility  of  application.  The  planned  segmented 
structure  of  the  program  also  makes  it  easy  to  insert  additional  cracked  and  conventional 
elements  as  desired. 

The  scheme  of  the  analysis  is  basically  divided  into  the  following  stages: 

•  Input  of  geometry,  material  data  and  restraints 

•  Input  of  applied  nodal  loads  and  imposed  nodal  displacements 
Calculation  of  element  stiffness  matrices 

•  Assembly  of  element  stiffness  into  a  system  stiffness 

•  Solution  of  simultaneous  equations  for  nodal  displacements 

®  Output  nodal  displacements 
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Calculation  and  output  of  forces  (or  stresses)  in  conventional  elements 

•  Calculation  and  output  of  stress- in  tensity  factors  and  strain-energy  release 
rates 

The  overall  logic  flow  of  the  computer  program  is  depicted  in  figure  5. 

The  computer  program  is  efficient  in  its  use  of  computing  time  and  core  storage.  Core 
use  can  easily  be  changed  to  fit  individual  problem  sizes  if  desired.  Full  advantage  has 
been  taken  of  symmetry  and  the  banded  nature  of  the  simultaneous  equations  to  be  solved, 
and  no  data  are  generated  when  not  required. 
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TABLE  II 


FUNCTION  OF  SUBROUTINES 


SUBROUTINE 

NAME 

FUNCTION 

INPUT 

Inputs  model  general  specifications 

CORIN 

Inputs  nodal  coordinates  and  restraints 

TRIN 

Inputs  triangular  elements  for  both  isotropic  and  anisotropic  cases 

SPRGIN 

Inputs  spring  elements 

AXIN 

Inputs  axial  elements 

CRKIN8 

Inputs  8-node  symmetric  cracked  elements  for  both  isotropic  and 
anisotropic  cases 

CRKNIO 

Inputs  10-node  unsymmetric  cracked  elements  for  both  isotropic  and 
anisotropic  cases. 

MATIN 

Inputs  material  properties 

DISPIN 

Inputs  imposed  displacements 

LOAD  IN 

Inputs  applied  loads 

CONNEC 

Generates  data  which  subroutine  BANDIT  requires 

BANDIT 

Renumbers  the  node  points  to  minimize  band  width 

BAND 

Determines  the  band-width 

KFORM 

Forms  a  complete  stiffness  matrix 

TRI 

Computes  triangular  element  stiffness  matrix 

SPRNG 

Computes  spring  element  stiffness  matrix 

AXI 

Computes  axial  element  stiffness  matrix 

CRAK8 

Switches  control  to  either  CRAK8I  or  CRAK8A 
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TABLE  II  (ConHnued) 
FUNCTION  OF  SUBROUTINES 


SUBROUTINE 

NAME 

FUNCTION 

CRAK8I 

Compul-es  stiffness  matrix  for  8-node  symmetric  isotropic  cracked 
element 

CRAK8A 

Computes  stiffness  matrix  for  8-node  symmetric  anisotropic  cracked 
element* 

CRAKIO 

Switches  control  to  either  CRKIOI  or  CRKIOA 

CRK10I 

Computes  stiffness  matrix  for  10-node  unsymmetric  isotropic  cracked 
element 

CRKIOA 

Computes  stiffness  matrix  for  10-node  unsymmetric  anisotropic  cracked 
element 

INVERT 

Inverts  a  matrix 

CPLXRT 

Computes  complex  roots  of  a  specified  4th  order  characteristic  polynomial 

CPWR 

Raises  complex  number  to  a  real  power 

THCOEF 

Transforms  thermal  coefficients  from  one  reference  axis  to  another 

THDISP 

Computes  thermal  displacements  for  both  symmetric  and  unsymmetric 
cracked  element 

CHOL 

Solves  simultaneous  equations 

OUTPUT 

Controls  output  data 

DISPOT 

Outputs  Displacements 

TRIOUT 

Calculates  triangular  elements  output 

SPROUT 

Calculates  spring  elements  output 

AXIOUT 

Calculates  axial  elements  output 
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TABLE  II  (Continued) 
FUNCTION  OF  SUBROUTINES 


SUBROUTINE 

NAME 

FUNCTION 

CRK08 

CRKOlO 

Calculates  Kj  for  the  8-node  symmetric  elements 

Calculates  Kj  and  Kjj  for  the  10-node  unsymmetric  elements 
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SAMPLE  PROBLEMS 


This  finife-element  computer  program  has  been  tested  extensively.  Both  symmetric  and 
unsymmetric  cracked  elements  have  performed  well  with  respect  to  accuracy  and  efficiency 
as  shown  by  sample  problems  presented  here.  These  results  were  achieved  with  a  relatively 
coarse  finite-element  grid.  With  refinements  in  the  grid,  even  more  accurate  results  would 
be  obtained.  The  symmetric  cracked  element  normally  gives  more  accurate  results  than  the 
unsymmetric  cracked  element,  which  is  understandable  in  view  of  the  fact  that  the  unsym¬ 
metric  cracked  element  has  fewer  degrees  of  freedom  that  it  can  bring  to  bear  on  the  first 
mode.  However,  the  unsymmetric  cracked  element  can  be  used  in  a  much  wider  class  of 
crack  problems  and  is  more  practical  for  industrial  applications. 

To  illustrate  the  capabilities  of  these  two  types  of  crack  elements  (symmetric  and  un¬ 
symmetric),  a  number  of  selected  sample  problems  are  included  in  this  report.  These  sample 
problems  were  chosen  to  demonstrate  (1)  the  accuracy  and  economy  of  the  elements,  and 
(2)  the  versatilities  of  the  elements  to  perform  analyses  for  structural  configurations  of  prac¬ 
tical  importance. 


Anisotropic  Case 


Center-Cracked  and  Double-Edge-Cracked  Orthotropic  Tension  Plates  -  A  center- 
cracked  and/or  a  double-edge-cracked  orthotropic  plate  (Figure  6)  subjected  to  uniform 
tension  as  solved  by  Snyder  and  Cruse  (ref.  17)  was  investigated.  The  plate  analyzed  is 
228.6  mm  long  and  76.2  mm  wide.  The  finite-element  models  used  for  both  the  center- 
cracked  plate  and  the  double-edge-cracked  plate  are  shown  in  figure  7.  The  half-crack 
lengths,  a,  analyzed  are  15.24  mm  and  22.86  mm  for  the  center- cracked  plate  and 
15.24  mm  only  for  the  double-edge-cracked  plate. 


Numerical  results  were  obtained  for  10  representative  graphite  fiber-reinforced  epoxy 
laminates:  (0)  ,  (iSO)  ,  Ct34)  ,  (t45)  ,  (+56)  ,  l+(^),  (90)  (90  /±45)  ,  (90  /±45)  , 

.  .  S  S  S  S  5  5  ^  ^  a 

(0/i45/90)  ,  Lamina  properties  used  were: 


=  144.795  GPa 


=  GPa 


E22=  11.722  GPa  ,  y^2=°-21 


Table  III  indicates  the  complex  roots  obtained  from  the  characteristic  equation  for  each 
laminate. 


Results  are  presented  in  Tables  IV,  V,  and  VI.  The  distribution  of  differences  from 
Snyder  and  Cruse's  results  are  summarized  in  table  VII.  Notice  that  the  results  generated 
using  8-node  symmetric  cracked  element  are  closer  to  Snyder  and  Cruse's  results  than  these 


33 


TABLE  III 


COMPLEX  ROOTS  OF  CHARACTERISTIC  EQUATION 


^2 

MATERIAL 

“1 

“2 

»2 

(O/i  45/90)^ 

0- 

0.9999 

0. 

1,0001 

{90^/-  45)^ 

0,2680 

0.7007 

-0,2680 

0.7007 

(%/-  45)^ 

0. 

0.5554 

0. 

0.8372 

(90)^ 

0. 

0,2704 

0. 

1,0522 

(-"45)5 

0.7763 

0.6304 

-0.7763 

0.6304 

(^30)3 

0.9648 

1 .0090 

-0.9648 

1 .0090 

(±60)3 

0.4951 

0,5177 

-0.4951 

0.5177  i 

(0>5 

0. 

0.9504 

0. 

3.6982 

(*  34)3 

0,9415 

0.8730 

-0.9415 

0.8730 

(*56)3 

0.5711 

0.5296 

-0.5711 

0.5296 

NOTE:  u.  -  O'.  +  jS, 


TABLE  IV 


STRESS-INTENSITY  FACTORS  OF 
CENTER-CRACKED  PLATES  (a  =  15.24  mm,  a/w  =  0.4) 


MATERIAL 

SNYDER  &  CRUSE 

K|  (MPa  -  /^) 

8-NODE 

CRACKED  ELEMENT 

lO-NODE 

CRACKED  ELEMENT 

•<1 

%  OFF 

%  OFF 

{0/-  45/90)^ 

1.670 

1.675 

40.32 

1.622 

-2.88 

1.679 

1.696 

4-1.01 

1.652 

-1.61 

{’04/i  45)5 

1.670 

1.702 

+1.92 

1.671 

40.06 

(90)3 

1.683 

1.669 

-0.83 

1.755 

+4.28 

(-45)3 

1.745 

1.634 

-6.36 

-7.34 

(-30)3 

1.712 

j 

1.643 

-4.03 

1.574 

-8.06 

(-60)3 

1.719 

1.683 

-2.09 

1.649 

-4.07 

(0)s 

1.645 

1.652 

40.43 

1.578 

-4.07 

TABLE  V 


STRESS-INTENSITY  FACTORS  OF 
CENTER-CRACKED  PLATES  (a  =22.86  mm,  a/w  =0.6) 


MATERIAL 

SNYDER  &  CRUSE 

8-NODE 

CRACKED  ELEMENT 

lO-NODE 

CRACKED  ELEMENT 

K|  (MPa  -  y  m  ) 

%  OFF 

%  OFF 

(0/-  45/90)^ 

2.396 

2.396 

40.01 

2.326 

-2.92 

(90/45)3 

2.422 

2.444 

40.91 

2.389 

-1.36 

(9O4/-  45)^ 

2.396 

2.446 

+2.09 

2.413 

40.71 

(90)3 

2.394 

2.419 

+1.04 

2.574 

+7.52 

(-"45)3 

2.599 

2.421 

-6.85 

2.393 

-7.93 

(-30)3 

2.518 

2.405 

-4.49 

2.302 

-8.58 

(-Yo)3 

2.533 

2.473 

-2.37 

2.430 

-4.07 

<°>S 

2.314 

2.336 

+0.95 

2.229 

-3.67 
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TABLE  VI 


STRESS- INTENSITY  FACTORS  OF 
DOUBLE-EDGE-CRACKED  PLATES  (a  =  15.24  mm,  a/w  =0.4) 


MATERIAL 

SNYDER  &  CRUSE 

8-NODE 

CRACKED  ELEMENT 

lO-NODE 

CRACKED  ELEMENT 

K|  (MPa  -  y  m  ) 

%  OFF 

%  OFF 

(0/-  45/90)^ 

1.722 

1.701 

-1.22 

1.649 

-4.24 

(90^^  45)3 

1.722 

1.729 

40.41 

1.685 

-2.15 

("V-  '*^■'3 

1.7C6 

1.725 

+1.11 

1.638 

-1.06 

(90)3 

1.709 

1.689 

-1.17 

1.752 

+2.52 

<-«)s 

1.700 

1.757 

+3.35 

1.766 

+3.88 

(^30)3 

1 .783. 

1.722 

-3.42 

1.697 

-4.82 

440)5 

1.717 

1.767 

+2.91 

1.736 

+1.11 

(0)3 

1.670 

1.660 

-0.60 

1.598 

-4.31 

434)3 

1.750 

1.737 

-0.74 

1.724 

-1.49 

456)3 

1.745 

1.765 

+1.15 

1.746 

40.06 
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obtained  from  the  10-node  unsymmetric  cracked  element.  Since  Snyder  and  Cruse  used  the 
boundary- integral  approach,  no  definite  conclusion  can  be  drawn  with  regard  to  the  abso¬ 
lute  accuracy  of  their  results.  Besides,  these  finite-element  solutions  used  a  large  integra¬ 
tion  step  size  (10  steps  between  nodes)  along  the  boundaries  for  the  numerical  integration  in 
forming  the  stiffness  matrix  of  the  cracked  element.  Accuracy  would  be  improved  by  re¬ 
ducing  the  integration  step  size  (increasing  the  number  of  steps  between  nodes). 

Longitudinally  Cracked  Orthotropic  Strip-  The  longitudinally  cracked  orthotropic  strip 
problem  as  shown  in  figure  8  was  solved  under  plane-stress  and  imposed  displacement  condi¬ 
tions.  An  analytical  solution  was  developed  based  upon  the  energy  equivalence  and  Sih  s 
results  of  ref.  18.  The  solution  is  given  below  as 


^\°11  °22“  °12^  / 


+  a 


66 


1/2 


(80) 


where  a-  =  Elements  of  the  material  compliance  matrix 

6  =  Imposed  displacement 

h  -  Height  of  the  plate 

Numerical  results  were  obtained  for  the  two  different  laminates  shown  in  table  VIII. 

In  one  case,  the  laminate  is  stiffer  in  the  direction  normal  to  the  crack.  In  the  other  case, 
the  laminate  is-stiffer  in  the  direction  of  the  crack.  Both  laminate  cases  were  solved  using 
a  single  cracked  element  and  the  two  different  model  representations  shown  in  figure  9. 

The  results  are  summarized  in  table  IX. 

45-Degree  Cracked  Finite  Orthotropic  Plate  -  This  problem  is  shown  in  figure  10.  The 
finite-element  model  (figure  11)  representing  this  plate  consists  of  130  nodes,  210  aniso¬ 
tropic  triangular  elements,  and  two  10-node  anisotropic  cracked  elements.  The  problem 
was  solved  under  uniform  remote  stress  and  plane-stress  conditions.  The  stress-intensity 
factors  obtained  are  given  below: 

K|  =  1.1376  X  10"^  MPa -/TiT  and  K||  =  1 .1784  x  10“^  MPa -/nT 

Both  Ki  and  Km  are  within  ±2  percent  (-1 .77%  and  +1 .75%,  respectively)  of  Sih's  solution, 
ref.  18,  (K|  and  K||  =  1.1581  x  10“ 3  MPa  -/m)  for  an  infinite  sheet,  and  they  satisfy 
NASA's  requirement  for  the  accuracy  of  the  program. 
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TABLE  IX 


SUMMARY  OF  RESULTS  FOR  A  LONGITUDINALLY 
CRACKED  ORTHOTROPIC  STRIP 


1  Cracked  Element  & 

1 17  Triangular  Elements 


12.474 


-0.02 


5.326 


-2.70 


Isofropic  Case 

Cracked  Tension  Plates  -  The  sample  problem,  shown  In  figure  12,  was  the  first 
problem  analyzed  with  the  symmetric  cracked  element.  It  exhibits  the  degree  of  accuracy 
which  has  been  consistently  achieved  in  numerous  subsequent  problems.  The  finite-element 
model  has  30  nodes,  35  triangular  elements,  and  one  8-node  symmetric  cracked  element. 
The  three  configurations  (the  single-edge  crack,  the  double-edge  crack,  and  the  center- 
crack)  were  all  individually  analyzed  with  this  one  model  for  an  a/w  ratio  of  1/3.  The 
model  grid,  which  is  quite  coarse,  results  in  the  single  -  edge-crack  model  having  57  dis¬ 
placement  degrees-of-freedom  (DOF),  while  the  double-edge-crack  and  center-crack 
models  have  51  DOF.  The  stress-intensity  factors  computed  using  these  configurations  are 
compared  with  ASTM  (American  Society  for  Testing  and  Materials)  values.  The  accuracy 
of  the  finite-element  predictions  are  impressive  (<1.5%  error)  for  all  three  cases.  Refine¬ 
ments  in  the  model  would  produce  steady  convergence  toward  ASTM  values. 

The  same  cracked  problems  were  solved  using  an  unsymmetric  cracked  element.  The 
finite-element  model  as  shown  in  figure  13  consists  of  54  nodes,  69  triangular  elements, 
and  one  10-node  unsymmetric  cracked  element.  The  results  are  not  as  good  as  those  ob¬ 
tained  using  the  8-node  symmetric  cracked  element,  because  the  8-node  element  has  more 
degrees  of  freedom  to  bear  on  the  first  node. 

An  eccentric  crack  problem  as  shown  in  figure  14  was  also  solved,  and  the  results 
compared  with  Isida's  solution. 

All  these  results  are  summarized  in  table  X  for  reference. 


Bi-Material  Cracked  Plate  -  This  problem,  as  shown  in  figure  15,  was  solved  in 
response  to  the  NASA‘s  Request  for  Proposal  (RFP  1 -.31 -4957),  dated  28  May  1974.  It  was 
required  that  the  stress-intensity  factors  of  this  isotropic,  generalized  plane-strain  problem 
must  be  within  jl2%  of  Erdogan's  results  for  an  infinite  plate  (ref.  19).  All  dimensions  and 
material  properties  are  given  in  figure  15. 


In  the  finite-element  model  shown  In  figure  16,  only  the  upper  half  of  the  plate  has 
been  modeled  due  to  its  symmetry  with  respect  to  the  x-axis.  It  contains  100  nodes,  148 
plane-strain  triangles,  and  2  eight-node  symmetric  cracked  elements. 


A  uniform  displacement  of  46.228  x  10  mm  was  imposed  on  the  top  boundary  nodes. 
This  was  arrived  at  as  follows.  For  plane  strain, 


H(1  - 


H(1  -  1^2  )  S2 


where  H  is  the  half-height  of  the  finite  sheet.  was  taken  as  unity. 
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SUMMARY  OF  TYPICAL  RESULTS  FOR  ISOTROPIC  TENSION  PANELS 


The  computed  stress-intensity  factors  for  the  two  crack  tips  are 


2.67496  X  10 


=  0.25067 X  10 


These  values  include  in  their  definition  a  factor  of  y  tr  .  To  compare  with  the  values  in 
the  Request  for  Proposal,  it  is  necessary  to  divide  byyiT~;  hence, 

=  1.50918  X  10“^  MPa  K2  =  0. 14143  x  lO"^  MPa 

These  values  are  quite  close  to  those  computed  by  Erdogan  and  Biricikoglu  in  ref.  19.  The 
percentage  of  error  is  -0.12%  and  -tO.75%,  respectively.  The  axial  stresses  and  displace¬ 
ments  along  the  crack  line  are  plotted  in  figures  17  and  18. 

Stiffened  Plate  -  The  rate  of  growth  of  a  crack  in  a  stiffened  plate  will  be  influenced 
by  the  presence  of  the  risers.  The  simple  stiffened  plate  configuration  shown  in  figure  19 
has  been  analyzed  to  determine  this  effect.  The  plate  is  loaded  by  a  stress  along  its  edge 
normal  to  the  crack.  Two  cracks  are  assumed  to  propagate  from  the  center  of  the  two  edges 
of  the  plate  toward  the  risers.  When  the  cracks  have  passed  underneath  the  risers,  they  are 
assumed  to  extend  up  the  riser  and  into  the  plate  at  an  equal  rate  until  the  riser  is  com¬ 
pletely  failed.  The  crack  is  then  assumed  to  continue  to  propagate  into  the  plate. 

Only  one  quarter  of  the  configuration  shown  in  figure  19  is  idealized,  the  remainder 
being  represented  by  symmetric  boundary  restraints.  Both  the  plate  and  the  riser  are 
idealized  using  triangular  membrane  elements.  The  8-node  crack  element  is  used  to  repre¬ 
sent  the  crack  tip. 

The  results  of  the  analysis  is  shown  in  figure  20  in  the  form  of  the  crack  tip  symmetric 
stress-intensity  factor  plotted  against  the  crack  length.  It  is  evident  that  the  riser  will 
have  a  considerable  effect  on  the  rate  of  crack  growth. 


Thermal  Stress- Intensity  Problem  -  Sih  (ref.  21)  has  given  the  following  plane-stress 
formulae  for  the  stress- intensity  factors  appropriate  to  a  crack  with  constant-temperature 
faces  (see  figure  21)  acting  as  sinks  for  the  steady,  radially  inward  flow  of  heat  Q  at 
infinity. 


E  O'  Q(1  +  v)J~^ 

A\J1^  k 


(82) 


In  (82),  a  is  the  coefficient  of  thermal  expansion,  and  k  is  the  thermal  conductivity. 
We  now  solve  the  associated  temperature-distribution  problem  to  obtain  thermal  input  in  a 
form  acceptable  to  the  computer  program. 


The  general  real  solution  for  the  temperature,  T,  in  a  steady  heat-conduction  problem 
in  two  dimensions  is 
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vanish  on  the  crack  faces;  i.e. 

f '  (x)  +  f  "  (x)  =  O  ,  <  a^  (83) 

In  (87)  and  (88),  prime  (')  denotes  differentiation  with  respect  to  the  parenthetical  variable, 
while  the  superscript  (+)  or  (-)  indicates  a  value  taken  in  the  upper  and  lower  half  planes, 
respectively. 

The  boundary  condition  at  infinity  is 


(89) 


in  which  r  and  0  (shown  in  figure  21)  are  defined  by 

ie 

z  =  re 

In  view  of  (89),  f'(z)  is  expected  to  vary  at  infinity  like  z 

We  now  make  the  usual  substitution  for  a  cut-condition  like  (88);  i.e.. 


(90) 


(91) 
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V  2  2  2  2  1  /2 

z  -a  is  fhat  branch  of  (z  -a)  '  varying  like  z  at  infinity.  Because  of 

the  anticipated  remote  behavior  of  f  (z),  we  expect  g(z)  to  tend  to  a  constant  value  at 
infinity. 


The  problem  case  in  terms  of  g(z)  is  then 


g  (x)  -  g  (x)  ~  0  ,  X  <  a  ; 

(92) 

1  i  m  /V 

g(z;  -  c, 
z  ®  ^  1 

(93) 

The  obvious  solufion  being 

II 

CO 

(94) 

which  leads  to 

f(z)  =  c^  In  (z  +'y^z^  -  a^  )  +  c^ 

(95) 

when  (91)  is  integrated. 

The  constant  C2  in  (95)  only  sets  the  temperature  reference,  which 
ified ,  Taking  C2  —  0  with  no  loss  in  general  applicability,  we  find 

is  thus  far  unspec- 

2c 

^  1 

Sr  r 

(96) 

for  large  z.  Upon  imposing  (89),  we  find 

_  Q 
'"l  4kjT 

(97) 

Thus,  from  (86),  (95),  and  (97), 

T(x,  y)  l^Xn  (z  a  )  +  ^n  (z  +'^  -  a  )J 

(98) 

and  the  constant  crack-face  temperature  is  given  by 

.p  _  Q  in  a 
c  2k7r 

(99) 

Referring  to  figure  22,  it  is  convenient  to  let 
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then  the  temperature  distribution  can  be  written  as 


T(x,  y) 


^n  <  (r  cos  0  r«  cos 


(100) 


+  (r  sin  9  +y  r2  sin  - 2 —  / 


Figure  23  shows  a  finil-e-elemeni-  representation  of  the  first  quadrant  of  the  problem  in 
figure  21 .  The  eight-node  symmetric  cracked  element  ABCD  represents  the  crack-tip 
neighborhood  in  the  finite-element  model.  Input  particulars  are  taken  to  be 


255.928’='K  and  a  =  114.3  mm 
k 


(101) 


leading  to  a  crack-face  temperature  from  (99)  of 


T  =255. 505®  K 
c 


(102) 


and  triangular-element  temperatures  (centroidal)  from  (98). 

The  computer  program  was  executed  twice  with  additional  materfal  properties: 

E  =  68.95  GPa  ,  v  =  0.2  and  a  =  1 .8  x  10  ^  m.m  ^  .K  ^  (103) 

In  the  first  execution,  the  temperature  of  the  cracked  element  was  taken  to  correspond 
uniformly  to  that  of  the  crack  faces  (102).  The  computed  stress- intensity  factor 

K|  =  4.9137  X  10’^  MPa (^04) 

is  14.95%  higher  than  the  analytical  value  of  4.2747  x  10  MPa  -x/nT.  This  is  not 
unexpected  in  view  of  the  fact  that  the  entire  cracked  element  is  taking  the  crack-face 
temperature,  a  condition  thermally  more  severe  than  the  analytical  temperature  distribu¬ 
tion  (98).  This  is  reflected  in  the  results  of  the  second  execution,  in  which  the  tempera¬ 
ture  of  the  cracked  element  was  taken  to  be  255.567®K,  which  corresponds  to  the  average 
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of  six  superimposed  triangles  (shown  dashed  in  figure  23).  The  stress- intensity  factor 
computed  for  the  second  execution  was 


K,  =  3.6325  X  10“'^  MPa 


(105) 


which  is  15%  lower  than  the  analytical  value.  It  is  clear  that  more  refinement  in  the 
crack-tip  neighborhood  would  lead  to  a  more  tolerable  discrepancy.  It  is  important  to 
point  out  that  such  refinement  will  not  be  necessary  in  more  routine  applications  where  the 
temperature  gradients  near  the  crack  tip  are  less  severe. 


48 


APPENDIX  -  ELEMENT  STIFFNESS  MATRICES 


Axial  Element- 


The  positive  sign  convention  for  the  axial  element  Is  illustrated  in  figure  A-1 .  For 
/entlon  k 

T 


this  convention  k  and  Pj.^  are  as  follows: 


1  1 
1  1 


{p4=-AEc«T 


Triangular  Membrane  Elcmeni' 

The  sign  convention  for  the  triangular  membrane  element  Is  illustrated  In  figure  A-2 
Note  that  six  displacements  or  degrees  of  freedom  are  present.  It  is  then  reasonable  to 
assume  the  following  displacement  field. 

u  =  a^+  a2X  +  a2y  " ''4  V  V 

The  displacements  of  the  element  can  then  be  expressed  in  terms  of  the  a.  coefficients  as 
follows: 


1 

"2 

1  ^2 

°2 

^3 

°3 

< 

< 

> 

''1 

1 

°4 

'^2 

1  X2 

°5 

1 

CO 

>s 

CO 

X 

_ 1 

or  in  matrix  notation 

{«}  =  H  W 

then 

H  ■  H"  W 
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where 


>'3  >'3 


[c]  = 


''32  ■'^3  "^2 


^2^3 


■>'3  >'3 


^32  ‘^3  ^2 


where  X22  ~  ^3  ~  ^2' 


The  strains  can  also  be  expressed  In  terms  of  the  displacement  coefficients: 


_  _ 

X  ox  2 


_  3v  _ 

€  -  r—  -  a, 

/  oy  6 


ou  ov 


xy  ay  ox 


-  a„  +  a 


0  10  0  0  0 


e^V=  00000  1 


0  0  10  10 


If  the  strains  and  stresses  are  constant  within  the  membrane  element,  the  strain  energy  can 
be  expressed  as 


S  =  ’  A^,  {c-r 


where  is  the  area  of  triangular  element.  Let  the  stress-strain  relations  for  the  anisotropic 
material  in  the  elemental  reference  system  be  written  as 
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=  [A] 

Then  the  strain  energy  can  be  expressed  in  terms  of  the  node  point  displacements  as  follows: 


‘■ivH’M'WWH 

And  the|^kjmatrix  for  the  triangular  membrane  element  is,  therefore, 

H- ['■']’[’]'  H  H  ['■'] 

If  the  multiplication  is  carried  out  and  if  the  elements  of  the  matrix  are  arranged  so  that  the 
convention  in  figure  A-3  applies  then  the  k-  elements  of  the  stiffness  matrix  are 


kn  y3^  -  2/2X22  A^ 3  + X22  A23 


— )  *^12  ^3  ^^13  ”  ^3  ^32  ^^^12 ^^33^  ^  ^32  "^23 

V  t  / 


— Vl3““^3^'^n  ■^y3^^3‘^'^32^  3  “  ^3  ^32  ^^33 


'^14  “^3  ^1 3  "^^3  ^3 2  "^^3^32 '^33"  ^3^32  "^23 


— -)^]5  ~^2^3^\3^^2^32^32 


_-?j  *^16  '^2  ^3  "^12  ^^2  "^32  "^23 


S2  ^3  ^33"  ^^3  ^32 '^23 ''’^32  ^2 

V  t  / 


*"23  "  "  >*3^  ^13  ""3  ^3  ^^33  ^3  '^32  ^12  “  ""3  "^32  ^^23 

f  / 
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^4  =  ->'3^  ^33  ^  >^3  ^^3  ^  ^32^  ^23  ”  ^3  ^32  ^22 

*^25  "  ^2  ^32  ^23  “  ^2  ^3  ^33 

*^26  "  ^2  ^32  ^^22  “  ^2  ^3  ^^23 

*^33  ^^3  "^11  "  ^^3  ^3 '^IS  ^^3  "^33 

*^34  “  ^3  '^13  ~  ^3  ^3  ^‘^12  "^^33^  "^^3  ^23 

*^35  "  "^2  ^3  ^13  “  "^2  ""3  ^^33 

*^36  "  ""2  ^3^12  ~  '^2  "^3  ^23 

•^44  ys  ^33  "  ^^3  >^3  ^^23  ^3  ^^22 

*^45  "“^2  ^3  "^23^  "^2  ^3^33 

'^46  =  "^2  ^3  ^22  ■'^2  ^3  ^23 

k  =  X  ^  A 

55  2  33 

k  =  X  ^  A 

56  2  ^23 

k  =  X  ^  A 
66  2  22 
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where  “  2  ^2  ^3 
""32  "  ''3  ■  ''2 

The  values  of  P^t  for  the  triangular  plate  element  can  be  computed  in  terms  of  the 
elemeni*  stiffness  matrix.  Note  that 

{p)  =[k]{p} + {pj 

|f|p}=0  and{p[  is  selected  as  the  thermal  displacement  for  free  expansion,  then  the  forces 
for  complete  restraint  are 

For  the  triangular  plate  element,  the  components 

u,=0 

“2  “  “*1 1  ”2 
''3  =  “n 


''3 '“22  ^>'3 

The  symbol  *  implies  lhat  the  thermal  expansiwi  coefficients  hove  been  rotated  into  the 
local  reference  axis  for  the  element. 


Spring  or  Fastener  Element 

Figure  A-4  illustrgtes  two  plates  connected  by  a  shear  pin,  and  the  direction  vectors 
for  the  system  of  forces  that  are  assumed  to  act  on  the  pin  are  illustrated  m  ^ 

the  spring  element,  it  is  necessary  to  define  four  node  points.  Node  points  (1)  and  (2) define 
the  element  itself.  In  most  structural  models  these  node  points  may  be  coincident  before  any 
strains  occur,  that  is,  (Zi  -  Z2)  =  0.  Node  points  (3)  and  (4)  ore  used  in  con, unction  with 
node  point  (1)  to  define  direction  vectors  for  spring  elements  K]  and  K2,  respectively.  K3 
is  assumed  to  be  a  spring  element  perpendicular  to  the  plane  containing  K,  and  K2.  In 
other  words,  the  system  of  spring  elements  is  that  illustrated  in  figure  A-6.  If  we  let  u,  v, 
w,  be  displacements,  the  element  forces  for  the  spring  system  are 
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"kI 

=  K,  (g, 

-“2> 

’’yl 

=  Kj  (V, 

-V2) 

'■zi 

CO 

11 

-  w^) 

P  „ 

=  -p 

x2 

x1 

P  „ 

=  -p 

y2 

yi 

p  „ 

=  -p 

z2 

zl 

Hence,  the  stiffness  matrix  is 

r  K,  0 


0 


0  0k 
0  0 
0  -K^  0 

0  0  -K 


0  K, 


0  0 


where  K, ,  K^,  and  are  linear  spring  stiffnesses  for  the  fastener  element. 

10- Node  Cracked  Element  Thermal  Effects 

figurel-/'  ’u"  in 

9  •  I  he  stress  strain  law  can  be  written  as 


°12  °16 
°12  °22  ''26 


°16  °26  °66j  I'^xy, 


‘^12^ 
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=  0 

''3 

-0!22T  a 

''l 

=  0 

"4  = 

-Oii^TA 

"2 

=  -a,^TA 

< 

It 

^2 

=  0 

CI!^^TA-CI!^2^  ^ 

"3 

=  -a^  ^  T  A  -  2^  ^ 

''5  = 

<3 

h** 

CN 

1 
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w^  =  q.,,TA 

Vg  A 

v^  =  0 

u^  =  -a^  ^TA  +a^2T^ 

^TA  +a^2'^^ 

'^9  “  ^22^ ^ 

u,o--a,,TA 

Ug-^ijTA 

''10°® 
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FIGURE  1 .  TRIANGULAR  ELEMENT  MATERIAL  AXIS 


FIGURE  2.  CRACK-TIP  NEIGHBORHOOD 


(a)  8-node  Symmetric  Cracked  Element 


(b)  lO-node  Unsymrnetric  Cracked  Element- 


FIGURE  3.  CRACKED  FINITE  ELEMENTS 
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FIGURE  4.  SEGMENT  STRUCTURE  FOR  COMPUTER  PROGRAM 


8- NODE  SYMMETRIC 
CRACKED  ELEMENT 


10- NODE  UNSYMMETRIC 
CRACKED  ELEMENT 


FIGURE  7.  CRACKED  ORTHOTROPIC  TENSION  PLATE 
FINITE-ELEMENT  MODELS 


CRACKED 

ELEMENT 


FINITE-ELEMENT  MODELS  FOR  A  LONGITUDINALLY  CRACKED  ORTHOTROPIC  STRIP 
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FIGURE  12.  SYMMETRIC  CRACKED  ELEMENT  MODEL  FOR 
SINGLE-EDGE,  DOUBLE-EDGE,  AND  CENTER 


t  f  HH  f 


UNIFORM  DISPLACEMENT,  EQ,  (81) 


YOUNG’S  MODULUS 
E^  ®68.950  GPa 

=  30.683  GPa 

POISSON'S  RATIO, 

=^0,30 

1^2  =0.35 

STRESS  INTENSITY  FACTOR 
(REFERENCE  19) 

=  1.375S^ 

=  2.768S2 

WHERE  S  IS  REMOTE  STRESS. 


figure  15.  BI-MATERIAL  CRACKED  PLATE 
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FIGURE  1/.  STRESS  IN  BI-MATERIAL  CRACKED  PLATE 


484  TRIANGULAR  ELEMENTS 
NUMBER  OF  DOF:  370 


FINITE  ELEMENT 
(NO  RISER) 


FIGURE  20,  EFFECT  OF  RISER  ON  STRESS  INTENSITY  IN  PLATE 


FIGURE  A-5.  DIRECTION  VECTORS  FOR  FASTENER  ELEMENT 
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FIGURE  A-7.  LOCAL  COORDINATE  SYSTEM  FOR 
TEN-NODE  CRACKED  ELEMENT 
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complete,  and  a  lasting  contribution  to  existing 
knowledge. 

TECHNICAL  NOTES:  Information  less  broad 
in  scope  but  nevertheless  of  importance  as  a 
contribution  to  existing  knowledge. 

TECHNICAL  MEMORANDUMS: 
Information  receiving  limited  distribution 
because  of  preliminary  data,  security  classifica¬ 
tion,  or  other  reasons.  Also  includes  conference 
proceedings  with  either  limited  or  unlimited 
distribution. 

CONTRACTOR  REPORTS:  Scientific  and 
technical  information  generated  under  a  NASA 
contract  or  grant  and  considered  an  important 
contribution  to  existing  knowledge. 


TECHNICAL  TRANSLATIONS:  Information 
published  in  a  foreign  language  considered 
to  merit  NASA  distribution  in  English. 

SPECIAL  PUBLICATIONS:  Information 
derived  from  or  of  value  to  NASA  activities. 
Publications  include  final  reports  of  major 
projects,  monographs,  data  compilations, 
handbooks,  sourcebooks,  and  special 
bibliographies. 

TECHNOLOGY  UTILIZATION 
PUBLICATIONS:  Information  on  technology 
used  by  NASA  that  may  be  of  particular 
interest  in  commercial  and  other  non-aerospace 
applications.  Publications  include  Tech  Briefs, 
Technology  Utilization  Reports  and 
Technology  Surveys. 


